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Aims. IRAS 19410+2336 is a young massive star forming region with an intense outflow activity. Outflows are frequently studied in 
the near-infrared (NIR) since the H2 emission in this wavelength range often traces the shocked molecular gas. However, the mech- 
anisms behind the H 2 emission detected in IRAS 19410+2336 have not been clarified yet. We present here spatially resolved NIR 
spectroscopy which allows us to verify whether the H2 emission originates from thermal emission in shock fronts or from fluores- 
cence excitation by non-ionizing UV photons. Moreover, NIR spectroscopy also offers the possibility of studying the characteristics 
5_( . of the putative driving source(s) of the H2 emission by the detection of photospheric and circumstellar spectral features, and of the 

"*T? ' environmental conditions (e.g. extinction). 

Methods. We obtained long-slit, intermediate-resolution, NIR spectra of IRAS 19410+2336 using LIRIS, the NIR im- 
ager/spectrographer mounted on the 4.2 m William Herschel Telescope. As a complement, we also obtained J, H and K s images 
with the Las Campanas 2.5 m Du Pont Telescope, and archival mid-infrared (MIR) Spitzer images at 3.6, 4.5, 5.8 and 8.0/um. 
Results. We confirm the shocked nature of the H 2 emission, with an excitation temperature of about 2000 K, based on the analysis 
of relevant H 2 line ratios, ortho-to-para ratios and excitation diagrams. We have also identified objects with very different properties 
and evolutionary stages in IRAS 19410+2336. The most massive source at millimeter wavelengths, mml, with a mass of a few tens 

0^ ■ of solar masses, has a bright NIR (and MIR) counterpart. This suggests that emission -probably coming through a cavity created by 

one of the outflows present in the region, or from the outflow cavity itself- is leaking at these wavelengths. The second most massive 

l/~) ' millimeter source, mm2, is only detected at A > 6pra, suggesting that it could be a high-mass protostar still in its main accretion 

r _| phase. The NIR spectra of some neighboring sources show CO first-overtone bandhead emission which is associated with neutral 

*— «, . material located in the inner regions of the circumstellar environment of YSOs. 
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1 . Introduction ciently high to overcome the radiation pressure and form mas- 
sive stars via continuous disk-accretion processes. Higher spa- 

H . IRAS 19410+2336 is a young massive star forming region which m reso i ution observations at 1 .3 mm, also with the Plateau de 

has been extensively studied in a series of papers during the Bure Interferomete] . show that these two dumps split up into a 

last few years. It is considered as a prime example of the com- eyen i arge r number of su b-structures, with at least 12 cores in 

plexity of high-mass star formation due to the intense outflow each one (Beuther & Schilke ^ 2004) These cores have masses 

T^L^° ^11 T ' ! u T ™ !t , C °, n „ tmUUm maP ° ranging between 1 .7 and 25 M Q , and visual extinctions of up to 

IRAS 19410+2336 obtained with the IRAM 30 m telescope at 10 00 magnitudes 

low spatial resolution reveals two massive adjacent star form- . . . , , , ....... 

, ul ,. ■ • ., ., j. .■ ■ At the very center of the southern clump and coinciding with 

ing clumps roughly aligned in a north-south direction and sep- , J . . , r^- — 7- , Jt -„ ,, 

arntpd hv -S3" (Rpnthprpt al 9007M Thp snnthprn rlnmn the m0St P romment of the mm Sources, |Beuther et al. (2002d) 

araieu uy ~oj) (oeuinerei ai., z\jvz.d). ine soumern ciump, . . TT _ . „.. <~, T t nn 1 • 1 • 1 1 

.., c , . OAr . ,. ■ ., ,, ,f detected H9O and Class 11 CH3OH masers, which are considered 

with a mass of about 840 M a , is more massive than the north- . , , • , • , , ,• ^ 

, .., f inn ha <~>u .• • r^r\i^ in to be closely associated with the earliest stages of massive star 

ern clump, with a mass of 190 M . Observations in CCX2-1) . ■L , , . .... *. 

... ,.„ , , .. , . , . , , . • formation. Only at this position is there evidence for a weak 

with 11 angular resolution show that each clump is associ- , ■', , r , , . _ _,,_. 

. j .., ?■ 1 11 ^a ■ .u t * j- 1 mJy unresoved (at angular resolution of 077) cm continuum 

ated with a bipolar molec ular outflow in the east-west direc- J ! „ ... — s — -r\ Bno9i 

tion dBeuther et aUl2002cl) . The southern clump is resolved into , " -. , " l ' L 7 — ™* 

at least four cores in the 2.6mm continuum using the Plateau BeutheLetalJ d2002a) reported the detection of hard X-ray 

de Bure Interfer ometer, where as the northern clump is resolved emission from a number of point sources in the vicinity of 

into two cores (iBeufher et all [2003). These authors disentan- IRAS 19410+2336. They conclude that this X-ray emission is 

gle at least seven (possibly even nine) bipolar outflows in this due to intermediate-mass pre-main-sequence Herbig Ae/Be stars 

region with accretion rates of the order of lO _4 M yr _1 , suffl- or their precursors. 

So far, only the kinematic distance to this object is known 

Send offprint requests to: N.L.Martin-Hernandez, leticia@iac.es and there exists a distance ambiguity, with a near distance of 
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2.1 kpc and a far di stance of 6.4 kpc ( Sridharan et all |2002). 
However, authors like Beuthe ret alj ([2002a) consider this source 
to be located at its near kinematic distance because the derived 
outflow parameters are unreasonably high if the f ar kinematic 
distance is assumed (see also Beut her et alll2002cl) . Other argu- 
ments that strengthen the selection of the near kinematic distance 
are the following . (1) Th e CS(2-1) profile for IRAS 19410+2336 
(Bronfma n et all [1996) is fairly narrow (with a FWHM of 
3.1 km/s) as compared with other sources that also have a com- 
plex outflow structure, for instance IRAS 1850 7+0121, with a 
linew idth of 6.8 km/s at a distance of 3.5 kpc (Shepherd et all 
2007). At a distance of 2.1 kpc, the 40" beam of these observa- 
tions is probably detecting only the mm clump and not the lobes 
observed in CO(2-l). At the far distance of 6.4 kpc, however, 
the red and blue lobes of the outflo w would have be e n det ected 
as wings in the CS(2-1) profile. (2) lBronfman et alj d 19961) ob- 
served four IRAS point sources with far-infrared colors typi- 
cal of ultracompact Hn regions, including IRAS 19410+2336, 
which have very similar (LSR) velocities. These sources are 
within a galactic extent of about 1.7°, or about 60 pc at a dis- 
tance of 2. 1 kpc. This is the typical size of a giant molecular 
cloud (GMC), and it is likely that these four sources belong to 
the same GMC. At the far distance of 6.4 kpc, on the contrary, 
there would have to be a series of GMCs with the same velocity, 
or a single GMC with an extent of 180 pc, which is unlikely. (3) 
IWatson et all d2003l) observed IRAS 19410+2336 in the HI 10a 
and H2CO lines in order to solve the distance ambiguity. No 
emission was detected in HI 10a, but they detected H2CO in ab- 
sorption with a (LSR) velocity of 22.6 km/s and a linewidth of 
4.7 km/s. No absorption is detected at larger velocities. 

We thus assume a distance of 2.1 kpc. At this distance, the 
infrared luminosity of IRAS 19410+2336 is 10 4 L . 

Outflows are frequently studied in the near-infrared (NIR), 
since the H2 emission in this range often traces the hot 
(< 2000 K) shocked molecular gas. In both low-luminosity and 
high-luminosity outflows, the H2 emission delineates single 
or multiple bow-shaped features associated with the leading 
edge of the CO emission and, in some cases, als o traces colli- 
mated jet s and cavity walls (e.g. see the review by Richer etall 
2000). iBeuther et alj d2003l) mapped the H2 emission around 
IRAS 19410+2336 with the 3.5 m telescope on Calar Alto and 
used this information (combined with high-spatial resolution CO 
observations) in their interpretation of the outflows occurring in 
this region. However, the excitation mechanism behind this H2 
emission has not yet been determined. With this aim in mind, 
we obtained long-slit, intermediate-resolution, NIR spectra of 
IRAS 19410+2336 using ORIS, the NIR imager/spectrographer 
mounted on the 4.2m William Herschel Telescope (WHT). 
These observations help us verify whether the H2 emission is 
produced by thermal emission in shock fronts or by fluorescence 
excitation by non-ionizing UV photons. These mechanisms can 
be distinguished since they preferentially populate different lev- 
els producing different H2 line ratios. NIR spectroscopy also of- 
fers the possibility of studying the characteristics of the putative 
driving source(s) of the H2 emission by the detection of photo- 
spheric and circumstellar spectral features. Moreover, it allows 
us to investigate environmental conditions (e.g. extinction) to- 
wards IRAS 19410+2336. 

As a complement, we also obtained J, H and K s images with 
the Las Campanas 2.5 m Du Pont Telescope. Combined with 
archival Spitzer data at 3.6, 4.5, 5.8 and 8.0 pm, obtained by 
the GLIMPSE Legacy Science Program (Benjamin et al., 2003), 
these images are helpful to detect optically obscured, deeply em- 
bedded young stars and protostellar candidates, excess emission 



from their circumstellar matter and, at earlier stages, from their 
inf ailing envelopes. 

This paper is structured as follows. Section [2] describes the 
observations and data reduction, and Sect. [3] presents the pho- 
tometric and spectroscopic results. Finally, Sect. [4] presents the 
discussion and main conclusions. 



2. Observations and data reduction 

2.1. Imaging 

J, H and K s images were taken with the Las Campanas 2.5 m 
Du P ont Telescope . This telescope is equiped with a NIR cam- 
era dPersson et all 1 1992b using a Nicmos-4 256 x 256 pixel 
HgCdTe array detector. The observations of IRAS 19410+2336 
were taken on 1996 May 26. Integration times of 5 seconds per 
frame were used. A dither pattern using 20-25" offsets were ap- 
plied resulting in an on-source integration time of 35 seconds. 

A color composite image of IRAS 19410+2336 is presented 
in the left panel of Fig. Q] Relevant stars are indicated. 

The measured FWHM is l'/2, l'/l and 1" in the J, H 
and K s band, respectively. We obtained the photometry of 
IRAS 19410+2336 using the standard IRAlQ DAOPHOT PSF- 
fitting routines on the final mosaicked frames. We constructed a 
semi-analytical PSF from 10 isolated stars in the field for each 
band. The PSF postage stamp radius was set to 3'.' 15 while the 
source fitting radius was 1". 2MASS zeropoints were estimated 
comparing the fluxes of 14 bright sources -selected from other 
fields observed during this same night- with their brightnesses 
in the 2MASS All-Sky Catalogue of Point Sources. We did not 
apply any color correction. A total of 47, 67 and 110 sources 
were extracted in the J, H and K s band, respectively. 

The electronic Table 1 lists the final J, H and K s magnitudes 
in the 2MASS photometric system for the stars in the field-of- 
view depicted in the left panel of Fig. Q] The final uncertain- 
ties are the result of propagating the errors coming from the in- 
strumental magnitudes and those from the linear fits used in the 
magnitude transformation. Limiting magnitudes are around 17.5 
mag in J, 17 mag in H and 16.5 mag in K s , 

Astrometry was performed by matching the positions of stars 
in common in our images and in the 2MASS survey catalogue. 
The astrometry accuracy is of the order of 077. 

Archival Spitzer images at 3.6, 4.5, 5.8 and 8.0 pm are also 
used. The contours of the 5.8 pm data are overplotted in the right 
panel of Fig. Q] Fig. [2] shows a more detailed comparison be- 
tween the Spitzer and the K s imaging data. Saturation affected 
the brightest components in the 4.5 and 8.0 pm band images, 
and so they are not shown. Saturation limits at 4.5 and 8.0 pm 
are, respectively, 200 and 740 mJy for a frame time of 2 sec- 
onds. Spitzer/IRAC magnitudes in the 3.6, 4.5, 5.8 and 8.0//m 
filters are also listed in Table 1 . They have been obt ained from 
the GLIMPSE Archive (e.g. iBeniamin et"aTl 12003b . available 
through the NAS A/IPAC Infrared Archivefl 3-<x detection limits 
are 15.5, 15.0, 13.0 and 13.0 mag in the 3.6, 4.5, 5.8 and 8.0 pm 
bands, respectively. 



1 IRAF is distributed by the National Optical Astronomy 
Observatory, which is operated by the Association of Universities 
for the Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation (http://iraf.noao.edu/) 

2 http://irsa.ipac.caltech.edu 
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Fig. 1. (Left) Las Campanas color composite image of IRAS 19410+2336 where J is blue, H is green and K s is red. Stars which are 
discussed in the text are labeled. (Right) In grey-scale we show the broad-band (K s filter) acquisition image of IRAS 19410+2336 
obtained with LIRIS. A dotted box outlines the field-of-view covered by the Las Campanas JHK S data. The black small circles 
indicate the positions of the H2 knots 1 and 2 detected in one of the two LIRIS slits (slit #1). The large dashed circles mark the 
positions and approximate sizes of the two massive star-forming clumps detected at 1.2 mm with the IRAM 30 m telescope. The 
white contours outline the Spitzer 5.8 fim emission (more de tailed MIR images o f the central source are shown in Fig. |2|. The red 
asterisks mark the positions of the mm sources identified bv lBeuther et al.l d 2003|). The oran ge filled triangle marks the position of 
the cm continuum sourc e associated with H 2O and CH3OH maser emission dBeuther et all 2002d). Blue squares mark the X-ray 
sources observed by Beuthe r et al.l (12002a). The green stars are 2MASS point sources of interest. The (0", 0") position in both 
images corresponds to R.A. (J2000.0)=19 h 43 m l I s . 18 and Dec (J2000.0)=23°44'0472 [see the electronic edition of the Journal for a 
color version of this figure]. 



2.2. Spectroscopy 

NIR spectra in the 1.9-2 A /zm wavelength range (using 
the mrk grism) were obtained during the nights of 2006 
August 31 and September 1 and 2 at the 4.2 m WHT using 
LIRIS, a long-slit intermediate-resolution infr ared spectrograph 
dAcosta Pulido et all 120031: iManchado et~afl H)04). LIRIS is 
equipped with a Rockwell Hawaii 1024 x 1024 HgCdTe array 
detector. The spatial scale is 0725 pixeL 1 , and the slit width used 
during the observations was 0775, allowing a spectral resolution 
R - A/AA ~ 1000. The slit (of length 4.'2) was oriented along 
RA.= 0° (slit #1) and 59° (slit #2) covering the brightest nebular 
spots (see Fig. [3), Weather conditions were relatively good, al- 
though with sparse cirrus. The infrared seeing during our obser- 
vations varied between 077 and 178 measured from the FWHM 
of the standard star spectra. 

Observations were performed following an ABBA 
telescope-nodding pattern, placing the source in two posi- 
tions along the slit. Offsets were 15" for slit #1 and 25" for 
slit #2. Individual exposures were taken with integration times 
of 300 seconds, adding up to a total on-source integration time 
of 80 minutes for slit #1, and 60 minutes for slit #2. In order 
to obtain the telluric correction and the flux calibration, the 
nearby AOV stars HIP 95487 and HIP 95560 were observed 
with the same configuration. The data were reduced following 
standard procedures for NIR spectroscopy, using IRAF and 
LIRISDR, the LIRIS data reduction package. After the flat-field 



correction, consecutive AB pairs of two-dimensional spectra 
were subtracted to remove the sky background. The resulting 
frames were then co-added to provide the final spectrum. Sky 
lines present in the science raw data were used to determine the 
wavelength calibration. The final uncertainty in the wavelength 
calibration is of the order of 1A. The resulting wavelength- 
calibrated spectra were divided by a normalized reference 
spectrum to remove the telluric contamination. This normalized 
reference spectrum was generated from the observed spectra 
of the A0 V standard stars, divided by a Vega model spectrum 
convolved with the actual spectral resolution (~ 23A). For 
this, we used the rou t ine xte llcor -general within the Spextool 
package dVacca et al.L [2003). This routine was also used for 
the flux calibration. We estimate an uncertainty in the flux 
calibration around 20% based on the comparison of the various 
spectra obtained for the standard stars. This large uncertainty 
probably reflect the fact that sparse cirrus were present during 
the time of observation. 



The contours of the (acquisition) broad-band image of 
IRAS 19410+2336 obtained in the K s filter (cut-off wavelengths 
are 1.990 and 2.310 yum) are shown in Fig. [3] The slits are 
overplotted. We identify several regions, 1 to 4 in slit #1 (with 
P.A.= 0°), and 5 to 9 in slit #2 (with PA.= 59°). These regions 
represent well differentiated parts of the H2 emission. Regions 3 
and 4 of slit #1 correspond to two H2 knots (1 and 2) found at 
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Fig. 2. Spitzer/IRAC 3.6 and 5.8 jum images of the central 
part of IRAS 19410+2336. White contou rs represent t he K^ 
emission. The mm sources identified by iBeuther et al.l (2003) 
are indicated by asterisks. The triangle marks the position of 
the cm continuum source. In these figures, north is up and 
east is to the left. The (0", 0") position corresponds to R.A. 
(J2000.0)=19 h 43 m im8 and Dec (J2000.0)=23°44'04'/2 [see 
the electronic edition of the Journal for a color version of this 
figure]. 



about 50" and 1 19" north of object nr56 (see the right panel of 
Fig.[T] where their positions are indicated). 

Peak positions and fluxes of the lines were measured by fit- 
ting a Gaussian profile. These are listed in Table [2] The large 
discrepancies between regions 1 and 5 (both include the source 
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Fig. 3. Contours of the broad-band (K s filter) acquisition image 
of IRAS 19410+2336 obtained with LIRIS. The positions of the 
slits and the distinct regions analysed in the text are marked. 
Regions 3 and 4, which correspond respectively to H2 knots 1 
and 2 (see the right panel of Fig. [TJ, are 50" and 119" north 
of nr56. The red stars indicate, from right to left, the positions 
of sources nr56, nr56b, nr56c and nr56d. In this figure, north 
is up and east is to the left. The (0", 0") position corresponds 
to R.A. (J2000.0)=19 h 43 m l 1 S T8 and Dec (J2000.0)=23°44'04'.'2 
[see the electronic edition of the Journal for a color version of 
this figure]. 



nr56) are due to the different size and orientation of the extrac- 
tion apertures. A line is defined as being detected if its peak in- 
tensity exceeds the rms noise level of the local continuum (typi- 
cally around 4.5 x 10~ 16 Wm -2 //m _1 ) by at least a factor of 3. 



3. Results 

3.1. The environment of IRAS 19410+2336 

As mentioned in the introduction, IRAS 19410+ 2336 is formed 
by two massive adjacent star-forming clumps (IBeuther et all 
l2002bh . The southern clump, with a mass of 840 M , is about 
four times more massive than the northern clump. The positions 
and approximate sizes of these unresolved clumps are marked in 
the right panel of Fig. Q] The positions of the millimeter sources 
observed bv lBeuther et al.l (1200 3) are also indicated in Fig.Q] to- 
gether with the locations of the unresolved cm continuum emis- 
sion detected within the southern clump (Sridharan et al., 2002), 



and the X-ray sources observed bv lBeuther et al] (l2002al) . 

The elongated emission feature seen in the K s image co- 



incide s with an H2 linear structure observed by IBeuther et al.l 
(120031) . This elongated H2 feature is strongly associated with the 
blue- shifted molecular em ission of a CO outflow -denominated 
C by IBeuther et al. (2003)- and outlines the same morphologi- 
cal structure, indicating that this is a highly collimated outflow. 
This outflow is probably driven by the millimeter source mml. 
This is the brightest mm source, with a mass of a few tens of 
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Fig. 4. fLeffl NIR color-co lor diagram of IRAS 19410+2336. The solid line represents the main sequence from lTokunagal C2000) 
and lMartins & Plea d2006l) . The dashed lines show the reddening band for main sequence colors. A reddening vector of 2 mag in 
the /f-band is represented by a thick arrow. Blue squares indicate sources associated with X-ray emission. Orange diamonds are 
objects which are associated with strong MIR emission. Finally, various 2MASS point sources of interest are plotted by green stars. 
We include 2MASS_4, the 2MASS point source associated with the H2 knot 2. (Right) NIR color-magnitude diagram at a distance 
of 2.1 kpc. Here, the dashed lines show our detection limits, and the dotted lines comprise the area where reddened O and early B 
stars at a distance of 2. 1 kpc might be found [see the electronic edition of the Journal for a color version of this figure] . 



solar masses (Beuther& Schilke, 2004), and is associated with 
the NIR object nr56. It also has a MIR counterpart (see Fig. |2j). 
Another object of interest is source nr71, which also has 
a MIR counterpart detected by Spitzer, although weaker, and 
has no counterpart at millimeter wavelengths. Finally, the mil- 
limeter source mm2, w hich is nearly as massive as mml 
(Beuther & Schilke, 2004), is associated with MIR emission but 
remains undetected at wavelengths shortward of 6 //m. 



3.2. Color-color and color-magnitude diagrams 

The electronic Table 1 contains JHK S and IRAC magnitudes 
for 116 stars detected in the ~ 75" x 75" field around 
IRAS 19410+2336. The two bright red objects seen in the left 
panel of Fig. Q] are stars nr56 and nr71, with K s magnitudes of 
11.21 and 11.02, respectively. Source nr56 is only detected in 
K s . 

The color-color (/ — H vs H — K s ) and color-magnitude (K s 
vs H - K s ) diagrams are shown in Fig. [4] Error bars are 1-cr. In 
each diagram we have plotted the theoretical main sequence at 
the assumed distanc e for the region , i.e. 2.1 kpc. Intrinsi c color s 
were taken from iTokunagal d2000t) and IMartins & Plea* (120061). 
Absolu te visual m agnitudes, M v , were taken from ICramerl 
(fl997h . iHouk et al.l (U997I) and IMartins & Plea 1 d2006b . We also 
include four sources (2MASSJ - 2MASS_4) outside the ~ 
75" x 75" field around IRAS 19410+2336 which have avail- 
able 2MASS photometry. These sources are plotted by green star 
symbols and identified in the right panel of Fig.Q] The large scat- 
ter in the color-color diagram could be due to, among other pos- 
sibilities, the different spatial resolution between 2MASS and 



the Las Campanas observations, and the presence of blobs of 
emission with no detected star. 

The dotted lines in the color-magnitude diagram comprise 
the area where reddened O and early B stars at a distance of 
2.1 kpc might be found. Such stars would ionize their own Hn 
regions. Among the sources detected in H and K s , several dis- 
play high extinction appearing as red or orange in the left panel 
of Fig.Q] These stars are nrl4, nrl7, nr55, nr64, nr72 and nr97, 
with Ak in the range 1^4 mag. Stars nrl4 an d nr55 are associate d 
with two X-ray point sources detected by Beuthe ret alJ (|2002a). 
Among this group of candidate O and early B type stars are 
2MASS_1, associated with a mm source (mm5), and 2MASS_3, 
associated with extended MIR emission. From all these sources, 
only nr55 and 2MASS_3 exhibit a NIR infrared excess as can be 
seen in the color-color diagram. 

Other bright K s sources are stars nr34, nr56, nr71 and 
2MASS_2. From these, only star nr71 shows a clear NIR ex- 
cess; stars nr34 and nr56 are not detected in J and H. 2MASS_2, 
on the other hand, is associated with extended MIR emission. 

There are not many stars with a NIR excess even though 
the X-ray data indicates the presence of Herbig Ae/Be stars in 
the field around IRAS 19410+2336. However, this type of stars 
might only start showing an infrared excess at wavelengths long- 
ward of about 3 jum. 

Fig.|5]shows a [3.6]-[4.5] vs [5.6]-[8.0] color-color diagram 
for the sources with identifica tions in all four IRA C bands, where 
we use the classification of iMegeafh et al.l J2004I) . The box en- 
closed by dashed lines shows the domain of Class II objects, 
i.e. stars with disks. The domain of embedded young stellar ob- 
jects (YSOs), i.e. of Class and I objects, is also indicated. This 
diagram is useful to identify infrared excess emission of the ob- 



N. L. Martin-Hernandez et al.: Spatially resolved NIR spectroscopy of IRAS 19410+2336 



served targets. Sources nrl7 and 2MASS_2 fall within the "disk 
domain", while others, such as 2MASS_1, nrl4, nr34 and nr64, 
have a larger infrared excess and fall within the domain of em- 
bedded YSOs. In particular, 2MASS.1, nrl4 and nr34 have X- 
ray counterparts. Finally, nr56, with [3.6]-[4.5]~1.85, could cer- 
tainly be an embedded YSO. By comparison, nr71, with [5.8]- 
[8.0]~1. 16, does not have a large infrared excess and might be 
already a Class II object. 

3.3. Spatially resolved NIR spectroscopy 

Figure [6]illustrates the spatial variation of the H2 1-0 S(l) tran- 
sition at 2. 121 /im. Moving along slit #1 and from south to north, 
H2 emission is detected at the position of objects nr56 and nr71 
(regions 1 and 2), and about 50" and 1 19" north of object nr56 
(regions 3 and 4, also denominated H2 knots 1 and 2, respec- 
tively). Across slit #2, the H2 emission extends over ~ 25". We 
have distinguished 5 different regions in slit #2, named 5 to 9. 
Moving from NE to SW, the H2 line emission is stronger in re- 
gion 9. After this region, the slit crosses over two bright blobs 
(which correspond to sources nr56d and nr56c, see below), and 
decreases smoothly afterwards. 

The left panel of Fig. [7] shows the nebular spectra of 
IRAS 19410+2336 extracted from regions 2, 3 and 4 of slit #1, 
and the integrated spectrum across slit #2. The most obvious fea- 
tures of these spectra are a large number of strong molecular hy- 
drogen transitions and a weak Bry emission line. Table[2]lists the 
identified lines, their central wavelengths and fluxes. The quoted 
uncertainties are from the line fitting procedure. 

We also extracted four spectra along slit #2 where an unre- 
solved underlying continuum was detected. In particular, at 0" 
(which coincides with object nr56), at 3" (which we will de- 
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Fig. 5. Spitzer/IRAC color-color diagram for sources with iden- 
tifications in all four IRAC bands. The thick vector on the left 
shows Ak - 2. We use the same symbols here as in Fig. |4] The 
dashed lines and labels for different types of protostars are based 
on the classification of Mege ath et al.l Q2004) [see the electronic 
edition of the Journal for a color version of this figure]. 
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Fig. 6. Spatial variation of the H2 1-0 S(l) transition at 2. 121 //m 
across slits #1 and #2. The origin corresponds to the (0", 0") 
position of Fig. [3] The different regions analysed in the text (1 to 
9) are identified. Regions 3 and 4 correspond to H2 knots 1 and 
2, respectively. The stars indicate the positions of sources nr56, 
nr56b, nr56c and nr56d. 



nominate hereafter nr56b), at 8"7 (nr56c) and at 13'.'7 (nr56d). 
Since no photospheric features are detected, origins for the con- 
tinua other than stellar must be considered. Moreover, a stel- 
lar origin is unlikely because of the relatively large extinc- 
tion towards IRAS 19410+2336 (see Sect. E3). These con- 
tinua can also be produced by (1) free-free and free-bound 
emission from ionized gas, (2) scattered light and (3) hot 
dust. The free-free continuum emission detected at radio wave- 
lengths contributes only marginally to the infrared continuum 
de tected in the A"-band. This was e stimated using equation 7 
of Roman-Lopes & Abraham (2004) and the free-free Gaunt 
factors of iHummerl d!988l) . In particular, for a flux density of 
1 mJy at 3.6 cm and a typical electron temperature of 10 4 K, 
the expected free-free flux density in the K-band is only about 
0.15 mJy, which translates into ~ 10~ 15 W irT 2 //irT 1 at 2.2 //m, 
less than about 1 % of the observed continua. Most likely, the 
continua are caused by hot dust as the spectra of e.g. nr56 and 
nr56b are very red. These red slopes cannot be explained by scat- 
tered light, whose effect becomes stronger towards shorter wave- 
lengths. 

The spectra of nr56, nr56b, nr56c and nr56d are illustrated 
in the right panel of Fig. [7] The positions of these sources within 
IRAS 19410+2336 are indicated in Fig. [3] Besides the nebular 
H2 lines and weak Bry, an important result is the detection of CO 
first-overtone bandhead emission lines in nr56b and nr56c. The 
spectrum of nr56d has a very faint continuum and is mostly dom- 
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using the values listed in Table [2] 



inated by H2 lines. These four spectra have been corrected for 
extinction using the values listed in Table |2j i.e. Ak - 4.0 mag 
for nr56 and Ak =1.6 mag for the other sources. The derivation 
of these values will be explained later in Sect. 13.41 It is interest- 
ing to note that the continuum slopes of sources nr56c and n56d 
are not very red after dereddening. The red continuum slopes of 
sources nr56 and nr56b, as mentioned previously, can be easily 
explained by emission of hot dust. 



3.4. H 2 excitation 

The NIR emission from H2 can be produced by either ther- 
mal emission in shock fronts or fluorescence excitation by non- 
ionizing UV photons in the Lyman- Werner band (912-1 108 A). 



These mechanisms can be distinguished since they preferentially 
populate different levels producing different line ratios. 

Non-thermal excitation mechanims readily excite the v - 2 
and higher vibrational states, whereby collisional transitions 
preferentially de-excite the v = 2 level in favor of the v = 1 level. 
Hence, the H2 1—0 S(l)/2— 1 S(l) ratio has on occasion been used 
to distinguish between these two processes, with a pure fluo- 
rescense v alue around 2 and a thermal value of order 10 (e.g. 
Davi es et al.L 120 03). However, for high density clouds («h ^ 10 5 
cm -3 ) and for an intense incident FUV radiation field with Go > 
10 4 (where Go is the incident FUV photon flux between 6 and 
13.6 eV measured in units of the local interstellar radiation field), 
collisional de-excitation of FUV-pumped molecules will ther- 
malize the lower energy levels and the 1-0 S(l)/2-l S(l) ratio 
will rapidly approach a thermal value (e.g. Burton et al., 1990; 
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Table 2. Line fluxes measured along slit #1 (P.A.= 0°) and slit #2 (P.A.=50°) in units of 10~ 18 W irT 2 . We also give significant line 
flux ratios and results such as the ortho-to-para ratios (0i, 02 and 03), the extinction (Ak) and the excitation temperature (r ex ). 









Slit #1: 


: P.A.= 0° 








Slit #2 


: P.A.= 59° 








Region 1 


Re; 


gion 2 


Region 3 


Region 4 


Region 5 


Region 6 


Region 7 


Region 8 


Region 9 


5+6+7+8+9 


1-0 S(3) a 


7.4+1.6 


10.1 


±2.2 


22.1+5.9 


10.8 + 3.7 


20.1 ± 1.6 


14.7 + 2.5 


75.3+4.7 


55.2 + 4.4 


169 + 42 


343 ± 22 


1-0 S(2) 


3.0+1.2 




<4 


4.1+0.9 


1.6 + 0.2 


5.7+ 1.0 


4.3 + 0.5 


16.0 + 0.5 


16.2 + 0.6 


38.7+1.6 


81.8 + 1.9 


3-2 S(5) 


<5 




<4 


< 1 


< 1 


<3 


<2 


0.7 ± 0.2 


0.6 ± 0.2 


1.6 + 0.3 


2.9 + 0.9 


2-1 S(3) 


<4 




<4 


1.2 + 0.3 


0.6 + 0.3 


1.5 + 0.7 


1.5+0.3 


4.4 ± 0.2 


3.7 + 0.2 


10.8 + 0.5 


22.0+1.1 


1-0S(1) 


10.8 + 1.2 


4.0 


>±0.9 


10.2 + 0.9 


4.4+1.4 


20.5 ± 1.0 


12.3 + 0.5 


47.0 ± 0.8 


47.0+1.5 


106 + 5 


237 + 5 


3-2 S(4) 


<4 




<4 


< 1 


< 1 


<3 


< 1 


< 1 


< 1 


0.9 ± 0.2 


1.3 + 0.9 


2-1 S(2) 


1.4 + 0.9 


2.7 


±0.7 


<2 


<2 


<5 


<3 


1.9 + 0.3 


1.1+0.2 


4.8 + 0.3 


7.8+1.7 


Bry 


<3 




<3 


< 1 


< 1 


1.6 + 0.7 


1.2 + 0.4 


1.6 + 0.4 


< 1 


< 1 


5.4+1.7 


3-2 S(3) 


<6 




<6 


< 1 


< 1 


<4 


<2 


< 1 


0.5 + 0.1 


2.5 + 0.2 


4.3 + 1.4 


1-0 S(0) 


2.5 + 1.6 




<6 


3.0 + 0.4 


1.6 + 0.7 


6.4+1.4 


3.6+1.0 


11.7 + 0.4 


12.0 + 0.5 


26.8+1.4 


60.6 ± 2.7 


2-1 S(l) 


<6 




<6 


1.8 + 0.5 


< 1 


2.4+1.0 


1.4 + 0.5 


4.4 ± 0.4 


3.8 + 0.2 


11.6 + 0.7 


23.8 + 1.8 


2-1 S(0) 


<7 




<6 


<2 


<2 


<5 


<3 


<3 


0.6 + 0.3 


2.9 + 0.6 


5.1+2.6 


3-2 S(l) 


<7 




< 11 


<2 


<2 


<8 


<6 


<4 


1.2 + 0.5 


<3 


7.5+4.1 


1-0 S(l) (b) 
2-1 S(l) 


> 1.8 




>0.7 


5.7+1.6 


>4.4 


8.5 ± 3.6 


8.8 + 3.1 


10.7+1.0 


12.4 + 0.8 


9.1+0.7 


10.0 + 0.8 


1-0 S(I) (b) 
3-2 S(3) 


> 1.8 




>0.7 


> 10.2 


>4.4 


>5.0 


>6.1 


>47 


94+ 19 


42 + 4 


55+ 18 


01 


3.2+1.3 




>0.6 


2.4 + 0.3 


2.2 ± 0.8 


2.8 ± 0.4 


2.6 + 0.4 


2.9 + 0.1 


2.8 + 0.1 


2.8 + 0.1 


2.9 + 0.1 


<h 


<4.1 




< 1.9 


> 1.1 




>0.4 


>0.6 


2.5 ± 0.3 


3.8 + 0.4 


2.6 + 0.1 


3.3 + 0.5 


<Pi 


















> 1.0 


2.5 + 0.5 


2.8+1.5 


A K 


4.5+0.4 






2.8 + 0.5 


3.4 + 0.6 


4.0 ± 0.2 


1.6 + 0.4 


1.5 + 0.5 


1.2 + 0.7 


1.3 + 1.1 


1.6 + 0.5 


Tex (K) 


3440 ± 265 






2300 + 115 


2295 ± 165 


2000 ± 40 


2270 ± 80 


2115+70 


2000 ± 70 


2390+ 110 


2135 + 60 



(a) This line falls in a region of poor atmospheric transmission; the indentification as 1-0 S(3) is therefore dubious and it is discarded from the 
analysis. (b) This line ratio has not been corrected for extinction. 



iHollenbach & Nattal [19951 iDavies et"afll2003l) . Conversely, the 
higher vibrational states are expected to retain a fluorescent pop- 
ulation. For instance, the 1-0 S(l)/3-2 S(3) ratio retains a value 
of approximately 8, and only for very high densities and a very 
intense FUV radiation field will this ratio approach a thermal 
value of 10-100. Therefore, the combination of the 1-0 S ( 1 )/2— 1 
S(l) and 1-0 S(l)/3-2 S(3) ratios should indicate the H 2 excita- 
tion mechanism. 

Table |2] lists the 1-0 S(l)/2-l S(l) and 1-0 S(l)/3-2 S(3) 
ratios across slit #1 (PA.= 0°) and slit #2 (PA.= 59°). Across 
slit #2, the 1-0 S ( 1 )/2— 1 S(l) ratio remains practically constant 
with an average value around 10. The 1-0 S(l)/3-2 S(3) ratio, on 
the other hand, has an average value of 55 and reaches a value as 
high as 100 in region 8. The combination of both ratios suggests 
the presence of shocks across this bulk of emission. This analysis 
is however not conclusiv e in the case of slit #1, even though the 
dense PDR models of Davi es et al.l Q2003) seem to match the 
ratios observed for region 3. 

An additional diagnostic is provided by the ortho-to-para ra- 
tio. At the moment of formation, molecular hydrogen can join 
one of two denominations: ortho (aligned nuclear spins) or para 
(opposed nuclear spins). The ortho-to-para ratio 0, based on 
the statistical weights of the nuclear spins, is 3:1. In outflow 
regions where shock excitation is the primary emission mech- 
anism, observations of vibrational emission lines typically re- 
veal or tho-to-para ratios for vibrationally excited states of about 
3 (e.g. Smith et al., 1997). However, as discussed at length by 
Sternberg & Neufeld ( 1999), the ortho-to-para ratio can take val- 
ues other than the intrinsic 3 if the excitation is caused by UV 
fluorescence. Indeed, ratios in the range 1.5-2.2 have been mea- 



sured inJ^DRs_(e.g. Ramsay etal., 1993; Chrv sostomou et aU 
ll993HShupeetalJ. ll998: Lumsden et ail feoOll) . 

From the 1-0 S(0), 1-0 S(l) and 1-0 S(2) line fluxes we 
measure the ortho-to -para ratio for the v = 1 states (0i) fol- 
lowing equation 9 of ISmith et al.l (1 19971) . In Table [2] we list the 
resulting 0i for each region. In general, these are consistent with 
a ratio of 3 suggesting a thermalization of the v = 1 states. 

Of course, <p\ does not tell the whole story. A value of the 
ortho-to-para ratio for the v = 2 states (02) is derived from the 
2-1 S(l), 2-1 S(2) and 2-1 S(3) transitions using equation 1 of 
iDavis et all (120031) . Although the values for 02 are less certain 
because the v = 2 transitions are in general weak, the results 
in Table [2] indicate that this ratio is in general close to 3. This 
suggests a thermalization of the v - 2 levels, although we find 
02 < 1 -9 in region 2, a value which is typical of PDRs. 

In a few cases we could derive a value for the v = 3 levels 
using the 3-2 S( 3), 3-2 S(4) and 3- 2 S(5) transitions. Using the 
methodology of ISmith et al.l d!997l) . the ratio from these transi- 
tions is given by: 



1/03 = 0.945(/ 4 // 5 ) - 467 (/ 4 //3) a533 



(1) 



The values we obtain for slit #2 are again close to 3 as can 
be seen in Table [2] 

A more sophisticated way of characterizing the H2 emission 
is by plotting the observed column density against the energy of 
the upper level. The measured intensity, /, of a given H2 line can 
be used to calculate the column density of the upper excitation 
level of the transition, which, for optically thin emission, is given 
by: 



N; 



4nAl 
Ajhc 



(2) 
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Fig. 8. H2 excitation diagrams for every region (except region 2, where only two useful lines are available). The lines represent the 
best fitting purely thermal single-temperature models. 



where Aj is the rest wavelength and A, is the Einstein A- 
coefficient taken from Wolnie wicz et al.l d!998l) . If collisional de- 
excitation is assumed to dominate, the H2 will be in LTE and the 
energy levels can be described by a Boltzmann distribution. The 
relative column densities of any two excitation levels can thus be 
expressed in terms of an excitation temperature r ex : 



N: 



Sj 



<E> - Ej) 



(3) 



where gj i s the degeneracy, E ,■ is the energy of the upper level 
taken from lDabrowskil (I 1 9841) and k is the Boltzmann's constant. 
The values of Aj, Ej, A /and gj for the lines detected in the spec- 
tra are given in Table |3] Clearly, plotting the logarithm of the 
ratio of the column densities divided by the appropriate statis- 
tical weights gives the excitation temperature from the slope of 
the data. 

In Fig. [8] we show plots of logNj/gj versus the energy. We 
have normalized the population distributions relative to that in- 
ferred from the 1-0 S(l) line and corrected for extincti on. A 
standard extinction law of the form A,\ oc A~ 1J (e.g. iMathisl 



Il990t iMartin & Whittei 1 19901) was used. For thermalized pop- 
ulations at a fixed gas temperature, the logN/g points should lie 
on a straight line in these excitation diagrams. If the value of 
A/i used to correct the line intensities from which the column 
density is calculated is wrong, then the scatter of the points will 
increase. The values of A^ that minimize x 1 an d the resulting 
excitation temperatures are reported in Table [2] Of course, these 
values of the extinction are estimated assuming a purely thermal 
single-temperature gas and thus, they must be taken with cau- 
tion. Temperatures range from 2000 to 3400 K, while the extinc- 
tion spreads from approximately A^ - 1 mag to At = 4.5 mag. 

For fluorescence emission, each vibrational level has Nj/gj 
lying along a separate "branch" and the rotational population 
within each level can be approximated by a thermal distribution. 
A curved line in this plot would therefore provide evidence for 
non-LTE processes. In the case of IRAS 19410+2336, a single 
temperature LTE gas fits rather well the v = 1 and v = 2 tran- 
sitions indicating that at least these two levels appear to be ther- 
malized at temperatures of the order of 2000 K. Unfortunately, 
this cannot help to discriminate between shocks and UV excita- 
tion. This is because in dense PDRs the lower H2 levels will be 
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thermalized as in shocks. However, the distinction can be made 
by looking at the v = 3 transitions, which will be predominantly 
excited by the PDR. 

When v = 3 transitions are detected (as is the case of slit 
#2), their column densities are only slightly underpredicted by 
the single component models. This suggests that a non-thermal 
mechanism might be partially responsible for the excitation of 
these high vibrational levels. However, the fact that 1-0 S( 1)/Bry 
is s» 1 strongly ar gues against a stellar UV excitation model (cf. 
Hatc h et aUl2005h . Moreover, the faintness of Bry emission and 
the measured excitation temperatures suggest that fast C-shocks 
might be present in IRAS 19410+2336. This type of shock has 
typical velocities of the order of 40 kms 1 and heats th e gas t o 
~ 2000 K producing strong H2 lines (e.g. iDavies et al.L |2000). 
The UV radiation field necessary to excite the v = 3 levels of H2 
can be produced in the dissociated apex of a bow-s haped shock 
dFernandes & Brand! 11995b iFernandes et all Il997l) . The excess 
emission observed for the lines with v = 3 would then arise from 
H2 fluorescence produced by Lya pumping of the low density 
(~ 10 3 - 10 4 cirT 3 ) pre-shocked gas, while the bulk H2 emission 
is excited in a bow C-shock. 



Table 3. We give the air wavelengths, upper energy levels, 
Eistein A coefficients and degeneracies of the H2 lines detected 
in IRAS 19410+2336. 



Transition 


*i 


E, 


A, 


«\ 




Gum) 


(K) 


(10- 7 s) 




1-0 S(3) 


1.9571 


8365 


4.21 


33 


1-0 S(2) 


2.0332 


7584 


3.98 


9 


3-2 S(5) 


2.0650 


20856 


4.53 


45 


2-1 S(3) 


2.0729 


13890 


5.77 


33 


1-0 S(l) 


2.1212 


6956 


3.47 


21 


3-2 S(4) 


2.1274 


19912 


5.25 


13 


2-1 S(2) 


2.1536 


13152 


5.60 


9 


3-2 S(3) 


2.2008 


19086 


5.63 


33 


1-0 S(0) 


2.2229 


6471 


2.53 


5 


2-1 S(l) 


2.2471 


12550 


4.98 


21 


2-1 S(0) 


2.3549 


12095 


3.68 


5 


3-2 S(l) 


2.3857 


17818 


5.15 


21 



3.5. The origin of the Bry emission 

Recombination radiation from hydrogen atoms is observed 
through the Bry transition (A = 2.166 /mi) in regions 5, 6 and 

7 (see Table |2|, and is present in the spectra of nr56, nr56b and 
nr56c ; 

Beut her et alj d2002dl) reported the detection of a weak 
(~ 1 mJy) radio continuum source towards nr56. Using equation 

8 of [Martin-H ernandez et al.l d2003l) and considering optically 
thin emission, a dust-free Hn region with such 3.6 cm emission 
will produce a Bry line flux of about 1 x 10~ 17 W m" 2 , Reddening 
this line flux (we use an extinction A% - 4, see Table |2), this 
translates into an observed Bry line flux of ~ 2.5 x 10~ 19 W m~ 2 , 
i.e. about 15% of the Bry line flux measured in region 5. Other 
origins might be necessary to explain the Bry line emission de- 
tected in this region. In fact, Bry emission can also arise in 
shocks. It has been detected, for instance, in the bases of jets 
from T Tauri stars and in the accretion flows around low mas s 
Class I sources (e.g. lFolha & EmersonlfeOOUlDavis et aUl2001h . 

The Bry emission detected in the spectra of nr56b (region 
6) and nr56c (region 7), however, probably have a different ori- 
gin. The detection of CO bandheads in these spectra suggests the 
presence of dense circumstellar material, possibly a disk, imply- 
ing that the Bry emission could have the same origin. It could be 
produced, for instance, in a disk wind caused by the interaction 
between the UV photons and the ionized upper layer of the disk. 



4. Discussion and conclusions 

The new NIR spectroscopic data presented here has provided, on 
the one hand, a complete analysis of the H2 emission and, on the 
other hand, direct evidence for the presence of massive YSOs 
and their circumstellar environment. 



4.1. Stellar content 

The NIR emission of IRAS 19410+2336 mainly concentrates on 
an elongated feature suffering an extinction Ay ~ 15 mag. This 
structure is ali gned with the blue-s hifted molecular emission of 
a CO outflow (Beut her et all 120031) . At the emission maximum 
of this almost linear structure, where the source nr56 is located, 



the extinction is even higher, with A v ~ 40 mag. This source 
nr56 is not only associated with the brightest MIR source, but 
also with the millimeter source mml, which has a mass of a 
few tens of solar masses (Be uther & Schilke , 2004). A spherical 
distribution of this cold material would lead to a very large visual 
extinction (Ay ~ 1000). No emission at near-infrared, or even at 
mid-infrared wavelengths, would thus be expected. However, the 
observation of NIR and MIR counterparts show that emission is 
leaking at these wavelengths. One explanation would be that the 
NIR and MIR emission is coming through a cavity created by 
one of the many outflows present in this region. The NIR and 
MIR emission could be emitted in the close environment of the 
forming stars hidden in mml and escape via the outflow cavity. 
Another possibility is that the IR emission is coming from the 
outflow cavity itself, i.e. from dust heated on the walls of the 
outflow cavity near the source (iDe Buizer & Minien 120051) . In 
this picture, the mm emission would be associated with the outer 
regions of the circumstellar matter. 

Source nr56 also exhibits a large [3.6]-[4.5] excess, which 
indicates that it might be an embedded Class 0/1 YSO. However, 
the detection of 3.6 cm continuum emission suggests in principle 
the presence of a recently ignited, and thus a more evolved star 
which has already formed an ultracompact H n region. Assuming 
for now that there is indeed an Ultracomp act H 1 1 regio n -it could 
probably be a Hypercompact H n region dKurtzLl2005l) - and that 
the radio emission is opticall y th i n, it w ould be consistent with 
an ionizing B2 V star dPanagial Il973l) . However, in the case 
of the high-mass protostar NGC 7538 IRS l. lDe Buizer &Minierl 
(2005) note that the centimeter continuum emission detected in 
this object could be the partially ionized emission from an out- 
flow, or even, could arise from a photoevaporated disk wind. 
Shocks associated with the supersonic infall of matter onto a 
protostellar disk are also sources of radio free-free emission. 
This last scenario can however be excluded since the models 
predi ct fluxes much lower than 1 mJy (Neufel d & Hollen bach. 
1996). Radio observations of the contiuum emission at other 
frequencies will be most helpful to determine the nature of 
this emission, i.e. whether it is optically thin or optically thick. 
Another clue to the nature of nr56 is the detection of H2O 
and Class 11 CH3OH masers (Beut her et all [2002d), which indi- 
cates a very early evolutionary stage. Both maser types, although 
they are on occasion spatially coincident with ultracompact Hn 
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Fig. 9. Spatial variation across slit #2 (P. A.= 59°) of the ortho-to- 
para ratio (<f>\), the extinction (A K ) and the excitation temperature 
(T ex )- In grey we show the average results from the integrated 
spectrum across this slit. Offsets are the same as in the top panel 
of Fig. © 



regions, are more often found in hot molecular cores and/or 
massive molecular outflows. In particular, H2O masers are be- 
lieved to be excited by collisional pumping with H2 molecules 
within sho cks associat ed with outflows a nd/or accretion (e.g. 
lElitzur et allll989UGaray & Lizanol [l999). Finally, its ZT-band 
spectrum only shows weak Bry and is practically featureless. It 
has a red continuum slope which is probably caused by hot dust. 

Sources nr56b and nr56c, detected towards the east of nr56, 
show weak Bry and CO first-overtone bandhead emission in 
their spectra. This is the first detection of CO bandheads in 
IRAS 19410+2336. These CO ban dheads are observed in about 
25% of all the massive YSOs fe.g.lChandleretall fl993l fl99l 
iBik & Thil l2004t iBlum et all |2004j) and are emitted in neutral 
material with temperatures between 2000 and 5000 K, and den- 
sities of about 10 10 cirT 3 . This hot neutral gas is probably located 
in the inner, dust-free regi ons of the circum stellar environment, 
relatively close to the star. IBik & Thi (2004) modeled the pro- 
files of the CO bandheads in 4 massive YSOs and suggested that 
they could be explained by a simple Keplerian disk, probably 
the remnants of massive accret ion disks. Source nr56b al so has 
a weak millimeter counterpart (Beuther & Schil kel 120041) . Both 
sources appear to be less embedded -and perhaps more evolved- 
YSOs than nr56. Source nr56d, on the other hand, with a spec- 
trum dominated by H2 lines, is presumably a knot or bow shock 
in the outflow. 

Another source of interest is nr71, located 9" north of nr56. It 
has not a large IR excess and might be already a Class II object. 
It has a featureless /T-band spectrum. The lack of Bry emis- 
sion could indicate that the circumstellar material is not (or very 
midly) ionized. 



Finally, the millimeter source mm2, located 4" north of nr56, 
is associated with MIR emission only visibl e at A > 6 ^m. This 
source, which is nearly as massive as mml (Beuther & Schilke, 
2004), could be a candidate for a high-mass protostar, the equiv- 
alent of Class O low mass YSOs which are still in the main ac- 
cretion phase. 

Other sources in the vicinity of IRAS 19410+2336 show in- 
frared excess at wavelengths longward of about 3 /vm. In particu- 
lar, sources 2MASS_1, nrl4, nr34 and nr64 have a large MIR ex- 
cess, suggesting that they might be embedded YSOs. 2MASS_1, 
nrl4 and nr34 are associated with X-ray emission. Interestingly, 
there are not many objects around IRAS 19410+2336 with a NIR 



4.2. The nature of the H 2 emission 

Based on the H2 line ratios, the measured ortho-to-para ra- 
tios and results of the excitation diagrams, we conclude that 
the H2 emission observed across slit #2 is clearly produced 
by shocks. Of course, the fact that the linear emission feature 
m apped by this slit is aligned with one of the outflows found 
bv lBeuther et al.l (120031) reinforces the shocked nature of this H2 
emission. 

Figure [9] shows the spatial variation of the excitation tem- 
perature and extinction across slit #2. The temperature is prac- 
tically constant with a value around 2100 K, increasing slightly 
up to 2400 K in region 9. The extinction Ak, on the other hand, 
is higher in region 5, where it reachs a value of 4.0, and keeps a 
lower and almost constant value around 1 .5 from regions 6 to 9. 

The conclusions for the regions observed across slit #1 are 
more ambiguous because of the lower number of H2 lines de- 
tected. It is thus difficult to discriminate between shocks and 
dense PDRs. The low values found for <p\ and <p2 in regions 2, 3 
and 4 seem to favor the case of dens e PDRs. However, the pres- 
ence of shocks cannot be excluded. iNeufeld et al.l (Il998l |2006) 
have measured ortho-to-para ratios in Herbig-Haro objects and 
found values significantly smaller than the equilibrium ratio of 
3. They conclude that these low ratios imply that the gas has 
not been warm long enough to reach equilibrium between the 
ortho and para states and that they are the legacy of an earlier 
stage in the thermal history of the gas when the gas reached 
equilibrium at a temperature < 90 K. They even place a con- 
servative upper limit of ~ 5000 yr on the period on which the 
emitting gas has been warm. This upper limit is based upon 
the expected timescale for the conversion from para- to ortho- 
H2. Hence, the low ortho-to-para ratios we measured in these 
three regions could indicate either the presence of a very recent 
shocked gas or of a dense PDR. 

4.3. Conclusions 

We obtained long slit, intermediate-resolution, NIR spectra 
of IRAS 19410+2336 along two preferential directions, i.e. 
PA.=0° and PA. =59°. As a complement, we also obtained J, H 
and K s images taken with the Las Campanas Du Pont Telescope, 
which we combined with archival Spitzer data. 

We have found objects with very different properties and 
evolutionary stages in IRAS 19410+2336. In particular, the NIR 
spectra show CO first-overtone bandhead emission which is as- 
sociated with neutral material located in the inner regions of 
the circumstellar environment of YSOs. The brightest source 
is nr56, which is probably the driving source of the strongest 
outflow. It is also the most massive source at millimeter wave- 
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lengths, with a mass of a few tens of solar masses, and its early 
evolutionary stage is confirmed by the presence of H2O and 
CH3OH masers. The fact that nr56 is detected in the near- and 
mid-infrared suggests that the assumption of a spherical distri- 
bution is probably not right. Such a spherical distribution of cold 
material would lead to a visual extinction of up to 1000 magni- 
tudes. Emission is likely leaking at these wavelengths, coming 
through a cavity created by one of the outflows present in the 
region, or from the outflow cavity itself. 

The second most massive millimeter source, mm2, is only 
detected at A > 6/zm, suggesting that it could be a high-mass 
protostar still in its main accretion phase. 

Another important conclusion is the confirmation of the 
shocked nature of the H2 emission, with an excitation temper- 
ature of about 2000 K, based on the analysis of relevant H2 line 
ratios, ortho-to-para ratios and excitation diagrams. H2 emission 
at 50" and 119" north of nr56 is also found, but we could not 
determine whether it originates in shocks or in a dense PDR. If 
it is produced in shocks, this emission could be associated with 
the northern millimeter clump, where at least three outflows are 
present (Beuther et al., 2003). 
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Table 1. Photometry towards IRAS 19410+2336. 



ID 


R.A. (J2000) 


Dec (J2000) 


./ 


H 


K s 


[3.6] 


[4.5] 


[5.8] 


[8.0] 


1 


19:43:08.92 


23:43:31.2 


14.24 + 0.11 


13.79 + 0.13 


13.58 + 0.11 


13.46 + 0.08 


13.59 + 0.15 






2 


19:43:08.37 


23:43:31.4 






15.61+0.13 


13.02 + 0.07 


12.53 + 0.13 


11.80 + 0.20 




3 


19:43:12.46 


23:43:32.4 


15.56 + 0.12 


13.23 + 0.13 


12.03 + 0.10 


11.27 + 0.07 


11.16 + 0.08 


10.56 + 0.17 




4 


19:43:09.52 


23:43:33.1 


14.51 +0.11 


14.26 + 0.14 


13.92 + 0.12 


13.87 + 0.07 


13.90 + 0.20 






5 


19:43:08.29 


23:43:33.5 






15.68 + 0.14 










6 


19:43:12.63 


23:43:34.1 




15.86 + 0.17 


13.95 + 0.12 


13.04 + 0.13 


12.71 +0.14 






7 


19:43:07.55 


23:43:35.4 


17.55 + 0.20 


16.40 + 0.19 


15.82 + 0.15 










8 


19:43:13.26 


23:43:35.6 




16.49 + 0.23 


15.49 + 0.18 










9 


19:43:08.48 


23:43:36.4 


17.75 + 0.21 


16.09 + 0.19 


15.68 + 0.15 










10 


19:43:13.40 


23:43:36.2 


16.66 + 0.18 


14.80 + 0.14 


13.78 + 0.11 


13.03 + 0.10 


12.93 + 0.14 






11 


19:43:13.63 


23:43:36.3 






15.33 + 0.16 










12 


19:43:11.00 


23:43:37.0 




16.17 + 0.21 


14.01 +0.11 


12.85 + 0.10 


12.45 + 0.11 






13 


19:43:12.20 


23:43:36.9 






15.51+0.15 










14 


19:43:11.49 


23:43:38.1 


14.58 + 0.11 


12.95 + 0.12 


11.69 + 0.09 


10.64 + 0.05 


9.95 ± 0.06 


9.15 ±0.07 


8.47 ± 0.05 


15 


19:43:07.81 


23:43:39.1 






15.65 + 0.18 


14.48 + 0.17 








16 


19:43:07.52 


23:43:39.5 


16.70 + 0.16 


15.41 +0.16 


15.24 + 0.13 










17 


19:43:10.72 


23:43:40.9 




13.17 + 0.13 


10.96 + 0.09 


9.66 ± 0.08 


9.40 ± 0.06 


9.00 ± 0.05 


8.95 ± 0.08 


18 


19:43:09.12 


23:43:41.4 






16.57 + 0.23 










19 


19:43:12.62 


23:43:41.8 




16.55 + 0.20 


15.20 + 0.13 


14.34 + 0.14 








20 


19:43:08.10 


23:43:42.3 






14.21 +0.12 


10.49 + 0.07 


9.45 ± 0.05 


8.54 ± 0.05 


7.79 ± 0.05 


21 


19:43:11.60 


23:43:42.8 


15.57 + 0.15 


14.12 + 0.14 


12.85 + 0.10 


11.79 + 0.07 


11.12 + 0.08 


10.15 + 0.21 


9.46 ± 0.08 


22 


19:43:10.53 


23:43:42.7 






15.61+0.23 










23 


19:43:09.26 


23:43:43.7 


17.53 + 0.35 


16.77 + 0.25 


15.77 + 0.17 










24 


19:43:12.41 


23:43:44.7 






15.43 + 0.19 










25 


19:43:11.46 


23:43:45.2 


11.40 + 0.08 


11.23 + 0.11 


11.09 + 0.09 


11.14 + 0.06 


10.87 + 0.13 


9.93 ± 0.20 


9.97 ±0.16 


26 


19:43:09.49 


23:43:46.0 


17.02 + 0.17 


15.34 + 0.15 


14.31 +0.12 


13.22 + 0.15 


13.05 + 0.30 






27 


19:43:12.20 


23:43:46.5 






16.79 + 0.31 










29 


19:43:11.86 


23:43:48.7 




16.17 + 0.18 


14.11 +0.12 


13.02 + 0.12 


13.02 + 0.34 


10.75 + 0.33 




30 


19:43:12.36 


23:43:47.8 


13.78 + 0.11 


12.99 + 0.12 


12.56 + 0.10 


12.06 + 0.06 


11.88 + 0.13 






31 


19:43:12.84 


23:43:51.0 






15.82 + 0.22 










33 


19:43:10.09 


23:43:53.6 


15.10 + 0.11 


14.62 + 0.14 


14.07 + 0.12 










34 


19:43:11.83 


23:43:54.2 






11.55 + 0.11 


9.74 ± 0.05 


9.12 ± 0.11 


8.52 ± 0.05 


7.91 ±0.05 


35 


19:43:08.80 


23:43:54.9 






16.71 +0.21 










36 


19:43:13.48 


23:43:56.0 






15.17 + 0.14 










37 


19:43:09.49 


23:43:55.8 


16.30 + 0.16 


15.54 + 0.16 


14.94 + 0.13 










38 


19:43:07.72 


23:43:55.6 






14.77 + 0.13 










40 


19:43:13.07 


23:43:56.1 






15.19 + 0.16 










41 


19:43:10.81 


23:43:56.2 


16.91 +0.17 


15.24 + 0.16 


13.67 + 0.11 










42 


19:43:10.12 


23:43:58.0 






14.28 + 0.17 










43 


19:43:13.26 


23:43:58.7 






16.21 +0.17 










44 


19:43:11.44 


23:43:59.1 


17.30 + 0.22 


15.62 + 0.16 


14.55 + 0.12 










45 


19:43:09.38 


23:44: 0.6 


11.52 + 0.08 


11.10 + 0.11 


10.87 + 0.08 


10.77 + 0.07 


10.76 + 0.32 






46 


19:43:11.84 


23:44: 1.4 


16.44 + 0.13 


15.10 + 0.15 


14.41 +0.12 
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Table 1. continued. 










ID 


R.A. (J2000) 


Dec (J2000) 


./ 


H 


K s 


[3.6] 


[4.5] 


[5.8] [8.0] 


47 


19:43:13.09 


23:44: 1.8 




15.63 + 0.17 


14.32 + 0.14 


13.06 ±0.17 


12.11 ±0.17 




50 


19:43:11.52 


23:44: 4.6 






14.66 + 0.23 








51 


19:43:10.54 


23:44: 4.3 






15.43 + 0.13 








52 


19:43:08.96 


23:44: 4.7 




15.67 + 0.16 


14.19 + 0.11 


13.48 ±0.10 


12.66 ± 0.33 




53 


19:43:08.28 


23:44: 4.9 


13.40 + 0.14 


12.38 + 0.12 


12.18 + 0.10 


12.14 ±0.10 


12.14 ±0.16 




54 


19:43:11.36 


23:44: 5.3 






13.46 + 0.15 








55 


19:43:10.78 


23:44: 5.3 


14.91 ± 0.03 


13.78 + 0.16 


12.21 ±0.10 








56 


19:43:11.19 


23:44: 4.3 






11.21 ±0.13 


5.32 ±0.17 


3.47 ± 0.22 


3.53 ±0.07 


57 


19:43:08.41 


23:44: 6.2 




17.02 + 0.53 


15.40 ±0.17 








58 


19:43:13.43 


23:44: 6.4 






15.25 ±0.13 








59 


19:43:11.46 


23:44: 6.6 






13.91 ±0.17 








60 


19:43:10.14 


23:44: 6.9 






15.23 ±0.13 








61 


19:43:11.74 


23:44: 8.3 






12.90 ±0.16 








62 


19:43:12.00 


23:44: 9.0 




15.37 + 0.19 


13.69 ±0.15 








63 


19:43:10.77 


23:44: 9.2 


15.77 + 0.13 


15.36 + 0.16 


14.95 ±0.13 








64 


19:43:08.99 


23:44: 9.4 


14.09 + 0.11 


12.92 + 0.12 


12.05 ±0.10 


10.21 ± 0.05 


8.96 ± 0.05 


7.74 ± 0.03 5.76 ± 0.03 


65 


19:43:08.47 


23:44: 9.5 






16.50 ± 0.20 








66 


19:43:13.03 


23:44: 9.2 


17.80 + 0.26 




17.08 ± 0.30 








67 


19:43:11.48 


23:44:10.6 






15.22 ±0.14 








68 


19:43:12.03 


23:44:11.3 






14.47 ±0.14 








69 


19:43:10.06 


23:44:11.4 






15.02 ±0.12 








71 


19:43:11.17 


23:44:12.1 




14.37 + 0.14 


11.02 ±0.09 


8.10 ±0.06 




5.37 ±0.07 4.21 ±0.11 


72 


19:43:13.27 


23:44:12.3 




13.62 + 0.13 


11.27 ±0.09 


9.63 ± 0.08 


9.32 ± 0.06 


8.83 ±0.28 


73 


19:43:12.26 


23:44:11.7 






13.56 ±0.16 








74 


19:43:11.60 


23:44:13.4 


17.63+0.29 


15.97 + 0.17 


15.31 ±0.13 








75 


19:43:08.90 


23:44:13.6 






16.48 ± 0.20 








76 


19:43:09.14 


23:44:14.4 


14.02 + 0.10 


13.41 ±0.13 


12.85 ±0.10 


12.08 ± 0.07 


11.60 ±0.08 


11.05 ±0.12 


77 


19:43:11.20 


23:44:15.7 






15.08 ±0.15 








78 


19:43:10.80 


23:44:15.4 






16.85 ± 0.27 








79 


19:43:13.27 


23:44:16.8 




16.57 + 0.21 


15.32 ±0.16 








80 


19:43:11.09 


23:44:16.8 






15.14 ±0.13 








81 


19:43:12.39 


23:44:17.2 




16.65 ± 0.29 


15.47 ±0.15 








82 


19:43:12.57 


23:44:17.9 




17.24 + 0.29 


16.58 ± 0.25 








83 


19:43:08.84 


23:44:18.5 




15.84 + 0.16 


14.70 ±0.12 


13.82 ±0.12 


13.53 ±0.21 




84 


19:43:10.76 


23:44:21.2 


14.30 + 0.10 


14.02 + 0.13 


13.75 ±0.11 


13.42 ±0.18 




11.10 ±0.24 


85 


19:43:12.44 


23:44:22.0 






16.74 ± 0.23 


13.55 ±0.20 






86 


19:43:11.12 


23:44:22.5 




16.39 + 0.30 


15.47 ±0.13 








87 


19:43:13.10 


23:44:23.6 






16.48 ± 0.34 








88 


19:43:13.22 


23:44:25.6 


18.65+0.32 


16.77 ± 0.22 


16.08 ±0.19 








89 


19:43:08.56 


23:44:25.9 




15.16 + 0.15 


14.15 ±0.11 


12.72 ± 0.09 






90 


19:43:13.64 


23:44:26.7 






15.91 ±0.25 








91 


19:43:13.45 


23:44:29.5 






15.40 ±0.14 


13.69 ±0.14 


13.24 ±0.13 




92 


19:43:11.07 


23:44:29.9 






16.12 ±0.18 









2 

& 
X 

o 
3 

:■-:, 
3 

p. 

o 



C/3 

T3 



< 

O 

a. 



70 



•8 



> 

in 



3' 

I 



Table 1. continued. 



ID 


R.A. (J2000) 


Dec (J2000) 


./ 


H 


K s 


[3.6] 


[4.5] 


[5.8] [8.0] 


93 


19:43:13.06 


23:44:30.2 






16.02 + 0.19 








94 


19:43:09.25 


23:44:32.9 






13.78 + 0.13 








95 


19:43:13.17 


23:44:33.8 






16.48 ± 0.21 








96 


19:43:10.81 


23:44:34.0 


15.94 + 0.13 


14.92 + 0.15 


14.42 + 0.12 


14.03 + 0.10 






97 


19:43:09.42 


23:44:34.4 


13.36 + 0.10 


12.53 + 0.12 


12.04 + 0.10 








98 


19:43:08.72 


23:44:36.0 






15.74 + 0.15 








99 


19:43:08.07 


23:44:36.0 


16.53 + 0.15 


15.35 + 0.15 


14.89 + 0.12 


14.39 + 0.12 


14.32 ± 0.27 




100 


19:43:08.94 


23:44:36.4 


14.90 + 0.12 


14.41 ±0.15 


14.35 + 0.12 








101 


19:43:07.54 


23:44:36.9 






15.70 + 0.13 


15.11 ±0.21 






102 


19:43:12.78 


23:44:36.7 






17.56 + 0.46 








103 


19:43:08.89 


23:44:37.4 






14.32 + 0.12 








104 


19:43:11.87 


23:44:37.6 






15.50 ± 0.20 


13.85 ±0.11 


13.62 ±0.18 


12.05 ±0.16 


105 


19:43:12.01 


23:44:40.6 






16.27 + 0.20 


14.07 ±0.11 


13.05 ±0.13 




106 


19:43:12.93 


23:44:41.0 






15.21 ±0.18 


13.92 ±0.07 


13.83 ±0.17 




107 


19:43:10.47 


23:44:43.2 




16.86 ± 0.23 


15.65 + 0.15 








108 


19:43:11.96 


23:44:44.0 




16.26 + 0.27 


15.82 + 0.15 








109 


19:43:07.64 


23:44:44.4 


14.53 + 0.11 


13.81+0.13 


13.57 + 0.11 


13.44 ±0.08 


13.33 ±0.16 




110 


19:43:10.04 


23:44:45.2 




16.38 + 0.18 


15.97 + 0.17 








111 


19:43:09.35 


23:44:46.2 




16.88 + 0.21 


15.54 + 0.17 








112 


19:43:12.61 


23:44:46.9 






17.07 + 0.26 


14.33 ±0.12 


13.72 ±0.17 




113 


19:43:07.70 


23:44:49.8 


17.03+0.18 


16.04 + 0.18 


15.63 + 0.13 


15.11 ±0.17 






114 


19:43:11.49 


23:44:50.4 






15.64 + 0.19 


13.80 ±0.12 


13.03 ±0.14 


11.97 ±0.19 10.74 ±0.13 


115 


19:43:12.38 


23:44:50.8 






15.27 + 0.19 


13.34 ± 0.09 


13.98 ±0.20 




116 


19:43:12.50 


23:44:52.2 




16.11 +0.22 


14.96 + 0.16 




13.65 ±0.14 
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